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CHAPTER  1 


BACKGROUND 

1.1  Introduction 

Large  portions  of  many  structures  consist  of  plates  and  panels,  which  are 
susceptible  to  flexural  vibrations.  The  vibrating  structures,  in  turn,  radiate  acoustic 
energy  from  their  surfaces  into  the  surrounding  medium.  Although  the  applications 
of  damping  mechanisms  would  reduce  the  vibration  level  and  hence  the  acoustic 
radiation,  the  study  of  acoustic  radiation  from  an  elastic  surface  still  plays  an 
important  role  and  constitutes  one  of  the  major  areas  of  concern  in  acoustics. 

Numerous  studies  have  been  devoted  to  acoustic  radiation  from  vibrating 
structures.  The  first  solution  to  the  radiated  power  from  an  elastic  plate  modeled  by 
the  classical  plate  theory  was  obtained  by  Skudrzyk  [1-2]  and  Heckl  [3-5]  for  a  time- 
harmonic  point  force,  and  by  Thompson  and  Rattaya  [6]  for  a  time-harmonic  point 
moment.  The  solution  for  the  acoustic  radiated  pressure  from  a  point-excited  plate 
using  the  classical  plate  theory  was  given  by  Gutin  [7],  Feit  [8]  and  by  Skudrzyk  [9], 
The  influence  of  fluid  loading  on  the  radiation  from  an  elastic  plate  was  investigated 
by  Maidanik  and  Kerwin  [10].  All  of  these  investigations  employed  the  classical  plate 
theory  which  fails  at  high  frequencies  where  the  phase  and  group  velocities  become 
infinite  as  the  frequency  or  the  wavenumber  becomes  unbound. 

To  improve  the  high-frequency  prediction  of  the  acoustic  radiation  from  an 
elastic  plate  excited  by  a  point  force.  Feit  [11]  employed  the  Timoshenko-Mindlin 
plate  theory  which  accounts  for  the  transverse  shear  and  rotary  inertia  of  the  plate. 
Unlike  the  classical  plate  theory,  the  Timoshenko-Mindlin  plate  theory  correctly 
predicts  the  high-frequency  flexural  phase  velocity  of  the  plate.  Feit.  however, 
confined  his  analysis  to  time-harmonic  frequencies  below  the  second  coincidence 


frequency.  Stuart  (12-13]  considered  the  same  problem  with  the  excitation  frequency 
extended  above  the  second  coincidence  frequency  and  thus  predicted  a  second  lobe  in 
the  radiation  pattern  due  to  the  thickness  shear  pole.  His  solution,  with  new  insights 
into  the  leaky  waves  emanating  from  the  plate,  is  more  accurate  when  the  angle  of 
observation  approaches  the  coincidence  angle. 

Since  the  leaky  wave  poles  may  approach  the  saddle  point,  the  first  order 
approximation,  obtained  from  the  standard  saddle  point  method  by  Feit  [11] 
and  Stuart  [12-13],  vanishes  at  the  grazing  angle  and  significantly  overpredicts 
the  radiated  nearfield  acoustic  pressure  near  the  coincidence  angle.  Seren  [14] 
recently  used  the  modified  saddle  point  method  to  accurately  predict  the  nearfield 
scattered  pressure  from  an  elastic  plate  with  line  force  and  line  moment  impedance 
discontinuities  due  to  an  incident  acoustic  plane  wave. 

In  this  thesis,  the  modified  saddle  point  method  is  used  to  predict  the  radiated 
acoustic  pressure,  which  is  expressed  in  terms  of  a  uniformly  asymptotic  series.  This 
asymptotic  series  solution  is  valid  in  the  farfield  and  the  nearfield,  and  is  especially 
accurate  at  the  grazing  angle. 

1.2  Basic  Approach 

The  structure  considered  in  this  thesis  is  an  infinite  fluid-loaded  plate  excited  by 
a  time-harmonic  force.  Although  it  is  known  that  the  infinite  plate  does  not  exist 
in  reality,  the  justification  for  this  assumption  lies  in  the  fact  that  the  structural 
damping  of  the  plate  is  sufficient  to  neglect  boundary-reflected  waves,  which  are 
assumed  to  be  of  relatively  small  amplitudes  as  compared  to  the  waves  propagating 
from  the  excitation  point. 

Throughout  this  study,  the  plate  is  assumed  to  be  homogeneous,  elastic  and 
isotropic.  Instead  of  using  the  classical  plate  theory,  the  Timoshenko-Mindlin  theory 
for  plates  (or  thick  plate  theory)  is  considered  which  includes  the  transverse  shear 
and  rotary  inertia  effects.  These  transverse  shear  and  rotary  inertia  corrections 


to  the  plate  equations  will  be  shown  to  provide  two  coincidence  effects  and  thus 
generate  two  peaks  in  the  acoustic  radiation  pattern  at  frequencies  above  the  second 
coincidence. 

The  wave  motion  in  the  acoustic  medium  is  derived  for  the  fluid  velocity 
potential.  At  the  plate-fluid  interface,  the  fundamental  boundary  condition  of 
continuity  of  normal  velocity  is  applied  to  the  coupled  plate-fluid  system.  Two  types 
of  excitation,  i.e..  line  and  point  force,  are  studied  in  this  thesis  and  the  radiated 
acoustic  pressure  is  obtained  from  three  different  methods,  which  originate  from  the 
same  integral  representation  along  the  steepest  descent  path. 

The  first  method  used  in  the  analysis  is  the  well-known  method  of  the  steepest 
descent  path,  also  known  as  the  saddle  point  method,  which  is  based  on  the 
assumption  that  all  the  singularities  of  the  integrand  are  located  far  away  from  the 
saddle  point  and  the  steepest  descent  path.  Consequently,  it  can  only  lead  to  a 
farfield  solution.  As  will  be  seen  in  the  following  chapters,  this  method  overpredicts 
the  radiated  pressure  in  the  nearfield,  especially  at  the  coincidence  angles.  In  order 
to  properly  treat  these  limitations,  the  saddle  point  method  is  modified. 

The  modified  steepest  descent  method  explicitly  isolates  all  the  singularities  of 
the  integrand  and  results  in  a  uniformly  asymptotic  series  solution  valid  not  only  in 
the  farfield.  but  in  the  nearfield.  The  grazing  field  can  be  accurately  predicted  by  this 
modified  saddle  point  method,  since  the  first  order  approximation  to  the  solution  by 
the  standard  steepest  descent  method  vanishes  at  the  surface  of  the  plate. 

The  third  method  involves  numerical  integration  along  the  steepest  descent 
path,  where  the  integrand  becomes  an  exponentially-decaying,  non-oscillating 
function.  Results  obtained  from  the  numerical  integration  are  compared  with  those 
predicted  from  the  previous  two  methods.  As  will  be  seen  later,  the  excellent 
agreement  between  the  uniformly  asymptotic  series  solution  and  the  numerical 
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The  acoustic  wave  equation  for  the  fluid  velocity  potential,  <£,  is  given  by 


V2$(f,t)  = 


1  d2$(f,t) 

C2  dt2  ’ 


where  c  is  the  phase  velocity  in  the  acoustic  medium. 

The  acoustic  pressure  is  related  to  the  velocity  potential  by 


p{r,  t)  =  po 


d*(r,t) 


where  po  is  the  fluid  mass  density. 

2.2  Line  Force  Excitation 

For  a  line-excited,  fluid-loaded  elastic  plate,  as  shown  in  Fig.  1.  the  total 
transverse  loading  can  be  expressed  as 


f(x,t)  =  F(t)S(x)  -  p(x,  z  =  0  +,t), 


■>> 


where  F(t)6(x)  represents  the  line  force  applied  at  x  =  0  and  the  rectangular 
coordinate  system  has  been  used. 

For  a  time-harmonic  dependent  line  force  excited  at  a  circular  frequency  ui. 

F(t)  =  FQe~ujt .  one  can  suppress  the  time  dependence  and  the  T-M  plate  equation 
of  motion  becomes 

[(V2  +  ^5cj2)(DV2  4-  rrxlu2)  -  mu2|w(i) 

D  (5 

=  (1  -  SV2  -  -SIu2)\F06(x)  -  p{x.  z  —  0+)|. 

Letting  $  =  the  acoustic  wave  equation  can  be  expressed  in  terms  of 


acoustic  pressure  as 


r  d 2  d 2 


where  kn  =  u n'c  is  the  acoustic  wavenumber  excited  at  frequency 
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Define  the  complex  Fourier  transform  pairs  as  follows: 


/-J-oo 

p(x,z)e~lkxdx ,  (?; 

-OO 

and 

1  Z"*"00 

p(z^)  =  —  J  p{k,z)e'kxdk.  (8; 

By  taking  the  complex  Fourier  transform  of  Eqs.  (5)  and  (6),  the  plate  equation 
of  motion  and  the  acoustic  wave  equation  become 


[(  — fc2  +  — 5oj2)(  —  Dk2  +  mluj 2)  —  mu2]w(k) 


=  (1  +  sk 2  -  ^S/w2)[Fo  -  P[k,z  =  0+)| 


The  solution  to  Eq.  (10)  can  be  obtained  as 


p{k,z)  =  Ae"Sk?>-kZz, 


where  A  is  the  wavenumber  (angular)  spectrum  function  for  the  plate.  In  order  to 
so've  for  the  unknown  constant  .4  in  the  transformed  /c-plane.  one  should  introduce 
the  continuity  equation  given  as 


—  (fc,2  =  O'")  =  p0u2u>{k), 
oz 


which  essentially  states  that  the  normal  velocity  at  the  plate-fluid  interface  must  be 
continuous. 

Applying  the  boundary  condition.  Eq.  (12).  to  Eq.  (11).  and  substituting  the 
transformed  displacement  tl'(/r)  into  Eq.  (9).  results  in 


w{k)  =  - -  x  k£  -  k2A, 

pnwJ-  V 


Oi 


n  =  —  is  the  normalized  frequency. 

uc 


C 


kc 


2  /m  2  / 12(1  -  u*)p  s/\2c2 
~  C  V  ~D  ~C  V  Eh*  “  cph 

CJc 

=  —  is  the  acoustic  wavenumber  at  wc. 
c 


is  the  classical  coincidence  frequency. 


cl  =  - - - —  is  the  phase  speed  of  the  plate  longitudinal  wave. 

p  (l-i/*)p 

For  the  Bernoulli-Euler  plate,  one  can  obtain  the  radiated  acoustic  pressure  by 
neglecting  the  transverse  shear  and  rotary  inertia  efFects.  or  simply  substituting  S  = 
I  =  0  in  Eq.  (17). 

2.3  Point  Force  Excitation 

For  a  point-excited,  fluid-loaded  elastic  plate,  shown  in  Fig.  1.  the  total 
transverse  loading  which  consists  the  external  force  and  the  acoustic  pressure  at  the 
plate  surface  can  be  expressed  as 


f{r,t)  =  ~  P(r,z  =  0+,t),  (18) 

where  F{t)  represents  the  point  force  applied  at  r  =  0  in  the  cylindrical 
coordinate  system. 

For  a  time  harmonic  point  force,  the  T-M  plate  equation  of  motion  and  the 

acoustic  wave  equation  become 

[(V2  4-  y^Su/2)( jDV2  4-  mlu2)  -  mu>2}w(r) 

D 

=  (1  -  SV2  -  ^5/w2)[Fo^  -  p(r,z  =  CT)!, 

U  <Jtt  r 

and 

[V2  -r  —  —  —  ko2',ptr,z)  =  0.  (20) 

az- 

The  Hankel  transform  is  applied  to  the  spatial  coordinate  r  and  the  transform 


pairs  are  defined  as 


G(k2,u2)  +  iDF(k2  ,u2) 


2  ,,,21  V  *»  ' 


p  l)W- 


G(k2,u2)  =  1  +  Sit2  - 


mu 


D 


■is, 


and 


F{k2  ,u2)  =  fc4  -  ^(J  +  S)k 2  +  )2/5  -  ^ 


Finally,  the  pressure  field  can  be  obtained  by  taking  the  inverse  Fourier 


transform  and  expressed  as 


Again,  to  eliminate  the  integration  along  the  branch  cuts,  the  transformation 
k  =  kr>sina.  r  =  Rsin4>,  and  2  =  Rcosd>  are  used  and  the  pressure  field  in  the 


spherical  coordinate  is  of  the  form 


p(R,d>) 


PqFq 

4nm.fl2 

(Ci  -f  Czk^sin2 a)H^  {knRsinosina)ke)Sina.cosa.elk"Rcoaac°3'i’ 
(F5  +  Fzstn2a)  +  icosa(sin4a  —  Fisin2a  +  F2) 


da, 


(29) 


where  F’s  are  shown  in  Eq.  (17)  and  the  constants  Ci  and  C 3  are  the  same  as  C0 
and  C 2  in  Eq.  (17).  respectively. 


CHAPTER  3 


SOLUTION  FOR  LINE  FORCE  EXCITATION 

3.1  Standard  Steepest  Descent  Method 

Consider  a  general  integral  of  the  form  I  —  fc^Qj  F[a)ep^a'1  da  with  a  large 
parameter  p.  where  /(a)  and  F(a)  are  arbitrary  analytic  functions  of  the  complex 
variable  a  and  C(a)  is  the  path  of  integration  in  the  complex  a-plane.  For  the 
integral  in  Eq.  (17).  one  has  p  =  k0r  and 

f(a)  =  icos(a  -  <t>)  =  fi[a)  +  1/2(0:),  (30) 

where  a  =  ar  +  fa,.  /i(a)  =  sm(ar  —  <j>)sfn/ia,  and  /2(a)  =  cos(ar  —  <p)coshax. 

The  path  of  integration  is  chosen  such  that  /]  has  a  maximum  value  at  some 
point  and  decreases  as  rapidly  as  possible  with  distance  from  this  point.  The 
location  of  this  point,  a  —  a  —  aT  +  fa,,  is  defined  by  the  root  of  the  equation 

~  I a=Q  —  -isin(a  -  o)  =  0.  (31) 

da 

The  solution  to  Eq.  (31)  is  a  —  <t>  =  nn.  for  n  =  0, 1,2-  •  •.  The  choice  of  n  =  0 
gives  the  saddle  point  in  the  range  —  7r/2  <  a  <  n/2  at  aT  =  0  and  a,  =  0.  and  the 
location  of  saddle  point  is  mathematically  equal  to  the  observation  angle. 

Since  /(a)  is  an  analytic  function  and  both  f\  and  / 2  are  harmonic  functions, 
i.e..  V2/j  =  0  and  V2/2  =  0.  f[a)  cannot  have  maxima  or  minima  in  the  entire 
a-plane.  Hence,  the  point  where  Eq.  (31)  is  satisfied  is  a  stationary  point.  The 
topography  near  a  would  be  a  surface  that  resembles  a  saddle,  i.e..  paths  originating 
from  the  saddle  point  a  could  either  descend  or  ascend. 


On  the  steepest  path  in  the  complex  plane,  the  imaginary  and  real  parts  of  the 
analytic  function  /(a).  /]  and  /2.  are  such  that  the  lines  of  most  rapid  change  in  f\ 


correspond  to  those  where  /o  is  constant.  Thus,  the  path  of  steepest  change,  either 
descent  or  ascent,  is  defined  by  the  equation 


/2(a)  =  /2(a)  =  constant. 


(32) 


or 

cos(ar  -  4>)cosf iq,  =  1.  (33) 

The  steepest  descent  path  (SDP).  shown  in  Fig.  2.  is  then  chosen  such  that 
/x(a)  —  /i(<5)  <  0  and  the  integrand  decays  exponentially  at  the  fastest  rate  as  the 
point  of  interest  moves  away  from  the  saddle  point  along  the  SDP.  The  orientations 
of  the  SDP  at  the  saddle  point  can  be  shown  to  be  — tt/4  and  3 7r / 4 .  It  is  convenient 
to  express  the  SDP  as  the  locus  of  the  points  given  by 

s2  =  f{<P)  -  /(a)  =  t  -  /(a),  -00  <  s  <  +00.  (34) 


Since  /(a)  always  contains  a  constant  phase  term  i  along  the  SDP.  the 
transformation  given  in  Eq.  (34)  makes  the  SDP  coincide  with  the  real  axis  of  the 
complex  5-plane  and  the  saddle  point  a  =  <p  correspond  to  s  =  0.  The  choice  of  s2 
instead  of  5  in  Eq.  (34)  insures  that  the  integrand  decays  exponentially  with  fcnrs2 
for  the  real  variable  5  ranging  from  -00  to  -‘-00. 

Using  the  notation 

F(a)£  =  G(5)’  (35) 

the  general  integral  then  becomes 


/  =  /  ep/(a)F(a)da  =  ep/i0) 
'C(a) 


£ 


e-',J'G(s)ds. 


(36) 


Again,  expanding  G'(s)  about  the  saddle  point  5  =  0. 


G(s)  -  (7(0)  ~G',0)5  -  -G"(0)52 


-...A.i-V.-1  -V-  »■*  -  •■•••  ,A«.wAaV  / 


and  the  integral  can  be  expressed  as 


/+oo  r  +oo 

e-^'ds  +  G'i 0)  /  e-p°2sds  +  •••}. 

-  oo  J  —  oo 


(38) 


For  large  p.  i.e.,  p  >>  1.  Eq.  (38)  can  be  approximated  by  keeping  only  the 
leading  term  as 

I  =  epf{<t,)G{  0)y/^>  (39) 

where  p  =  kor  and  /(<£)  =  i.  Thus,  the  first  order  approximation  of  the  integral  is 
obtained  as 

/  ... 

(40) 


where 


m  = 


Is  DP  *=  \/ 


(Cq  +  Cik^sin2  <i>)cos<p 


(F$  +  F3sin2<j>)  +  icos<j>(sin4<p  —  Fisin24>  +  F2) 

The  total  radiated  pressure  field  has  the  form 

p(r, <£)  =  +  27ri  V' res  (  between  C(a)  and  SDP  )},  (42) 

2nmu‘i 


and 


where 


res(a;)  = 


iV(a)elfc',rc0S^Q~^^ 
2F3sinacosa  +  iD(a) 


N(a)  =  cosa(CQ  +  C^k^sin2 a) 


(43) 


and 


D{a )  =  —  sm5a  -t-  4sm3Qcos2a  +  F1sm3Q  —  2F\sinacos2a  -  F^sma. 

Note  that  there  are  10  poles  in  the  complex  a-plane.  but  only  those  located 
between  the  original  integration  path  and  the  SDP  would  contribute  to  the  radiated 
acoustic  pressure.  Figure  2  shows  the  original  integration  path  and  the  SDP.  The 


lL 


lL 


typical  locations  and  the  frequency  loci  of  the  a-poles  are  given  in  Figures  3  and  4. 
respectively. 

3.2  Modified  Steepest  Descent  Method 

In  the  standard  steepest  descent  method,  the  integral  for  the  acoustic  pressure 
is  approximated  by  taking  only  the  leading  term  of  the  Taylor's  series  expansion 
about  the  saddle  point.  However,  it  is  not  an  acceptable  approximation  if  the  poles 
approach  the  saddle  point,  because  the  Taylor's  series  expansion  converges  only  in 
the  region  with  a  radius  of  convergence  equal  to  the  distance  from  the  saddle  point 
to  the  nearest  singularity.  Therefore,  the  standard  steepest  descent  method  must 
be  modified  so  that  all  the  singularities  of  the  integrand  are  explicitly  isolated.  The 
asymptotic  series  solution  obtained  from  this  modified  method  is  valid  in  the  whole 
field,  while  the  saddle  point  contribution  leads  to  the  farfield  solution  only. 

The  total  pressure  is  the  sum  of  the  integration  along  the  SDP  and  the  pole 
contributions.  Thus. 

,  ,  _  PqFq  f  (Cn  4-  C2kosin2a)cosaelk,,rco3la~0^ 

P  1  27rmf22  J sdp  (Fs  ■+•  Fzsin2a)  —  icosQ(jm4a  -  Fisin2a  —  Fz)  ° 


residues  of  poles  crossed  by  SDP  )}. 

(* 

As  mentioned  earlier,  the  steepest  descent  path  is  the  locus  of  points  given  by 
X*  =  f(d>)  -  f(a )  =  i  —  f(a ),  where  -  cc  <  \-  <  oc. 


Thus. 


X*  =  -i(cosa'  -  1)  =  2tstn2  — . 


where  a'  =  a  -  0  and  \2  is  real-valued  along  the  SDP.  The  following  identities  are 
obtained  from  Eq.  (45): 


v 


co.sq'  =  1  -4-  ix*,  sina'  —  \\  \2  -  2z.  and  da'  —  — -  ■  ■  d\  (46) 


Now.  the  integral  along  the  SDP  can  be  written  as: 


Is D p(ri  d)  =  v  2e 


/-t-oo  , 
1 


N(x)  e~k"rx~ 


where 


D(x)  =  D0  +  iD i(x)  +  ixVx2  ~  2isind>D2{x), 

Do  =  F3  -r  F$, 

D\  =  g5t'x10  +  ?4X8  +  93*X6  +  <72  X 4  +  g  1 1 X 2  +  9o, 

<75  =  cos3  0  4-  5  sin*  4>cos4>  —  10stn2<pcos3<p, 

94  =  5 cos3  4>  4-  25 sinA <t>cos<p  —  oOsin2  <frcos3  <t>, 

93  =  —  10cos5<j!>  —  40  sin4<t>cos0  +  90  sin2d>cos30  4-  3(Fi  —  2  )sin20cos0 

—  {F 1  —  2  )cos3<p, 

92  =  —  lOcos50  —  20sin4 cj>cos4>  4-  70stn2<£cos3  <£  4-  9(Fi  —  2)sin20cos0 

-  3(F\  -  2)cos30  4-  i(  —  Fzcos2^  4-  Fzsin2#), 

91  =  ocos3d>  -  20 sin24>cos3<t>  -  6(Fi  -  2)sin2  0cos0 

4-  3(Fi  -  2)cos3(£(l  4-  F2  -  F\)cos0  -  i2(- F3COS2 0  4-  Fjstn2#), 
90  =  cos3 0  4-  (Fi  —  2)cos3<£  4-  (1  4-  F2  —  F\  )cos0  +  1  (  Fjcos2^), 

D2(x)  =  PaX*  +  P3iX°  +  P2X4  +  P 1  * X=  -  Po, 


P4  =  -sm4o  -  5co54o  -+-  10sin2©cos2c>. 

P3  =  4  sin  4  0  —  20cos4o  -  40sin2ocos2o. 

Pi  =  4stn4o  -  30cos4o  -  5Ostn*0cos20  4-  3(Fi  -  2)cos20  -  (F 1  -  2)sm20. 

Pi  =  -20cos4o  -1-  20sm2ocos2o  -  6(Fj  -  2)cos2o  -  2(F|  -  2)si'n2o  -  :2 F^coso, 
pn  =  -5 cos’o  -  3(Fi  -  2)sin2o  -  ( 1  -  F;  -  Fi )  -  !2F?,co.«;o. 

Ar(X)  =  ^ 3 1  \  J  *  n2\4  ~  ri  1 1 X *  ~  rin  -  \  \  x2  -  2i{m2\A  -  m]!\2  -  m,,). 
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n3  =  C2&o(—  3sm‘  qcoso  4-  cos34>), 
n2  =  3C2k2(  —  3sin2ocoso  —  3  cos3d>), 
ri[  =  Cqcosq  ~  C  2k2  [7  sin2  ocosq  —  2  cos3q), 
no  =  C2k2sin2  qcosq  +  Cncoso , 
m2  =  C2/c2(  — 3cos2<!>s:n<z>  -4-  sin3o), 
m  1  =  C2/c2  (6co5*  tpsino  —  2sm3o), 

mo  =  —  Co-sm®  +  C2/c2(2cos2<2>.sm<i>  —  stn3<£). 

In  Eq.  (47),  the  poles  can  be  isolated  by  expanding  the  integrand  into 
partial  fractions.  The  x-poles  can  be  found  directly  from  the  a-poles  through  the 
transformation  given  in  Eq.  (45).  However,  it  is  not  easy  to  expand  the  integrand, 
since  the  denominator  contains  the  square-root  term.  \/x 2  -  2:.  To  eliminate 
the  square-root  term  and  to  expand  the  denominator  into  a  polynomial  of  \.  it  is 
necessary  to  multiply  both  the  numerator  and  the  denominator  by  D' (\)  given  as 

D'(x)  =  D0  4-  tDi(x)  ~  »X  V  X2  -  2isinoD2(x)-  (■ 


Eq.  (47)  then  becomes 
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J-°°  D(x)D’{x)y  1  4- 


where 


D(X)D‘(x)  =  £mx2)". 

n=  1 

X{x)D‘(x)  =  -4(x,o)  ~  \v  X2  ~  2iB(\.o), 


-4  ( X ,  o)  =  As(\,o)  -  1 1 
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perform  the  integration  in  a  closed  form.  To  treat  these  integrals,  let  b  =  -i\n  a 
H  -  \/  knr.  then  the  integrals  become 

roo  e~^x' 

P?(».\)=  - =dX,  ( 

Jo  (x2-62)Vi-T 


^M.x)  =  f 

Jo 
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The  square-root  term  can  be  expanded  into  Taylor's  series  about  the  saddle  point 
X  =  0.  and  expressed  as 


,/K  ill 

V  2 


00  -1-3  /1  m  • 
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Thus.  P*(ii,b)  and  P£(/j.,b)  can  be  written  as 


It  can  be  shown  that,  by  the  method  of  integration  by  parts,  the  second  integral 


in  Eq.  (55)  and  (56)  is  equai  to 


x7,62~  dX~{-\)kb-k  j 

7—1 


■2u-;k-> 


As  can  be  seen  from  Eqs.  (55-57).  only  one  type  of  integral  is  left  and  of  the 


form  [16] 
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0  X* 


d\  =  —erfc(bn)  e *  b  for  SR{6}  >  0. 
O'4  lb 
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Finally,  the  integral  for  the  radiated  acoustic  pressure  along  the  steepest  descent 
path  becomes 

^5Dp(r- <Z>)  =v  2e" fc,’r  *' 

M---(fc-i) 

k\ 
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(59) 


and  the  total  pressure  field  can  be  expressed  as 

pnFr, 


p(r,  o)  —  — - — %{/s  opfr-  ©)  —  2m  V'  resi  between  C(a)  and  SDP  j},  (60) 

!  9rrrrjD  -  * 
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where  the  residue  contributions  are  given  by  Eq.  (43). 

Instead  of  calculating  the  complementary  error  function  in  Eq.  (59).  the 
computer  algorithm  takes  advantage  of  expanding  the  w-function  into  an  asymptotic 
series.  The  following  two  equations  [16]  have  been  used: 


w|:  =  e  er 


fc(-iz)  =  ^ 


( i-P 


— •  n  -  -  i 

n  =  0M2 


.  for  all  c. 


(61 


and 


w(iz\  —  cz  erfc\z I  = 


.  for  J?{c}  >  0. 


(62 


For  an  argument  of  the  u>-function  greater  than  1.0,  continued  fraction 
representation.  Eq.  (62).  is  used  while  the  Taylor's  series  expansion.  Eq.  (61).  is 
used  for  the  argument  less  than  1.0.  Thus,  if  the  argument  of  the  w-function  is 
of  the  value  near  1.0.  then  the  series  converges  very  slowly.  Consequently,  many 
more  terms  have  to  be  considered  in  order  to  obtain  a  better  solution.  It  should  be 
noted  that,  in  general,  this  uniformly  asymptotic  solution  obtained  from  the  modified 
steepest  descent  method  runs  much  less  computer  CPU  time  than  that  for  the 
numerical  integration  to  be  discussed  later. 

For  the  farfield  solution.  Eq.  (59)  can  be  simplified  to  a  more  compact 
form.  The  asymptotic  expansion  for  the  tu-function  with  a  large  argument  can  be 
expressed  as 


;(2)  ss  -L-:l  +  (- 


l-3,5--(2m  —  1),  . 

1]  — (i^r - 11  or 


Thus.  Eq.  (59)  becomes 


Isop(r.o)  =  -  v  2e'^fc"r-  * 1 


10  oc 
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>■  sL*  ■  2(m*l|  L*. 

n  =  1  rn  =  0  Xn  k=  1 


1  •  3-5-  --(2k  -  1)  k<  ~  zfer.T (m  +  +) 

k\  l)  x2(m-fc  +  l)l(A.orjm+l  I' 

(64) 


For  large  /tor.  the  solution  can  be  obtained  by  keeping  only  the  leading  term 
(m  =  0)  and  expressed  as 


/  \  PoFn  2.T  ,lk  ,,r-~ 

PsDp(r-o)  =  kTP 


y  an  _  an 

(x2  -  Xl)\x=</ 


It  can  be  shown  that  Eq.  (65)  is  essentially  an  alternate  version  of  the  pressure  field 
contributed  at  the  saddle  point.  \  =  0. 

3.3  Numerical  Integration 

The  numerical  integration  method  is  discussed  in  this  section  and  the  results 
are  compared  with  those  obtained  from  the  previous  two  methods.  It  can  be 


UW  JV  JL  k  A.  .  A-'.  *U 


J  J 


1 


seen  from  Eq.  (15)  that  the  numerical  integration  in  the  complex  /C-plane  is  very 
complicated  due  to  the  branch  cuts.  However,  one  can  alternately  perform  the 
numerical  integration  along  the  steepest  descent  path  without  loss  of  accuracy. 
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As  discussed  in  the  modified  steepest  descent  method,  with  the  transformation 
X2  =  -t[co5(a  —  0)  —  1 ) .  the  steepest  descent  path  is  then  transformed  to  the 
real  axis  of  the  complex  x-plane.  The  integrand  is  given  in  Eq.  (49).  Since  the  SDP 
has  the  characteristics  that  the  phases  are  stationary,  the  integrand  now  becomes  an 
exponentially-decaying,  non-oscillating  function. 

It  should  be  noted  that  the  computer  CPU  time  for  the  integration  along  the 
steepest  descent  path  is  much  less  than  that  needed  in  the  complex  K- plane  which 
requires  integrations  on  the  branch  cuts. 


* 

% 

s 

s 

h 


23 


CHAPTER  4 


SOLUTION  FOR  POINT  FORCE  EXCITATION 


4.1  Standard  Steepest  Descent  Method 

The  integral,  as  shown  in  Eq.  (29),  can  be  evaluated  asymptotically  by  the 
saddle  point  integration.  The  asymptotic  representation  of  Hankel  function  is  given 


r/O)/  \  ^  /  jLpt(2-’r/4)|i  4 _ L_ 

m0  ^  i1  +  8t2 


128(12) 


For  the  farfield  solution,  taking  only  the  leading  in  Eq.  (66)  and  substituting  into  Eq. 
(29)  results  in  the  following  equation  which  is  ready  for  standard  steepest  descent 
method: 


p(R,d>\  =  ■■  —  / - £ - e-‘*V4 

Airmfl2  V  nkoRsiruj) 

f  (Ci  -r  C3k^sin2a)kosinacosaelk"Rcoa(a~o) 

*  /  - - - - - - - ; - ==  d(Xn 

JC{a)  [(^5  +  Fjsin2a)  -r  icosa{sin4a  -  Fism2a  —  F2)j\/5in<z> 

where  R2  =  r2  A-  z2.  As  was  done  in  the  line  force  excitation,  the  saddle  point  is 
found  mathematically  equal  to  the  observation  angle.  The  first  order  approximation 
to  the  integral  in  Eq.  (67)  can  be  obtained  as 

Isdp*  \^Rel(k''R'*)F{0)'  (4C 


where 


.  ICi  —  Csk^sin2 o)knsinocos0 

F(o)  =  — - — - 7 - 7 - p - - - - (68) 

j(rs  —  f‘3sm‘o)  —  icoso(sin*o  -  r\sin-o  —  r2) \\  stno 

Thus,  the  total  pressure  field  becomes 

._  ,  P^Fq  —2  ie'k"R  (Ci  —  Cqkr.xiTT  o)k,,sinocoso 

P\R,0)  ~  - - 3 - ; - - 

4rrmfi‘  knRsino  [F 3  -  Fjsin-o)  -  icoso{sinAo  -  Fisin-o  -  F2) 

-1-  2~t  residues  of  poles  between  5  DP  and  C(a)  }, 
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re5(QJ )  =  r - 7^—  \a  =  a,- 

zr^sinacosa  +  iU{  a) 


=  kosinacosa(Ci  4-  C^k^sin2  a)HQ  (koRsinosina) , 


and  D(a)  is  shown  in  Eq.  (43). 


4.2  Modified  Steepest  Descent  Method 


The  integral  along  the  steepest  descent  path  in  complex  a-plane  can  be 


expressed  as 


(Ci  4-  CzkrtSin2  a.)HQX\koRsinosina)kQsinacosaelk"Rco3ac°3,t’ 


(F$  4-  Fssin2a)  4-  icosa(sin4a  —  Fisin2a  4-  F%) 


To  proceed  with  the  solution,  it  is  convenient  to  use  the  integral  representation  of 


the  Hankel  function  given  as 


1%)  =  U»  [ 
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o  v  4 is  -  y 


Thus. 


IsDP  =  J 


f  (Ci  4-  C3k2sm2a)kosinacosaelk"lic03{'a 

sdp  (T5  4-  Fusin'2  a)  4-  icosa(sm4  a  —  F\sin2a  -r  F2) 


y  4iu  —  y: 


:dyda. 


where  u  =  koRsinosina.  Substituting  the  transformation  and  the  identities.  Eqs. 


(34)  and  (46).  into  Eq.  (73)  results  in 
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AT(x,<i>)  =  -2ix2{  —  C  ikosinocostp 


4-  C3k3\3(sin0cos3 d>  —  sin3 <j>cos<p) (1  +  t'x2)2  ~  sinocos3 0x2  (x2  —  2 1 ) ] } 


•  [(^3  +  Fs)  +  t(95*X  +  ?4X  +  <73tX  +  <7  2X  +  <7ilX  +  <?o)j 

—  2x2 sin<i){C iko(cos2 4>  -  st'n2<b)(l  +  t’x2) 

4-  Czk3\{Zsin2 4>cos2 4>  -  sin4d>)(l  4-  t'x2)3 
4-  (cos4g!)  -  3sin2<j>cos24>)(l  +  t'x2)x2(x2  -  2 1 ) ] } 

•  [P4X8  +  P3iX6  +  P2X4  +  Pi»X2  +  Pol, 


At{x,<p)  =  \Cik0{cos2<p  -  sin2<p)(  1  4-  t'x2) 


4-  C3kQ(3sin24)cos24>  —  stn4^)(l  4-  t'x2)" 


4-  C3k3(cos4<t>  -  3sin2(f>cos2<p)(l  4-  t’x2)x2(x2  ~  2i) ] 


•  \{F3  +  Fs)  +  t(gsi'x  +  g4X  +  <73  tx  +  <?2X  +  9i»X  +  9o)| 

4-  C3/Co[3(sm<ji)cos30  -  5tn30co5oi)(l  4-  t'x2)2  -  sin0cos3 <j>x2 (x2  —  2t)|} 
•  [(F3  4-  Fs)  +  i{gsixl°  +  ?4X8  +  93tX6  +  ?2X4  +  <7i*X2  +  9o)| 

+  {Ci Ar05tn(^co5<^[(l  4-  t'x2)2  -  X2(x2  ~  2*)] 


4-  C3k3{3sin0cos3d)  —  3sin30cos0)(l  4-  t’x  )  X*(x  —  2t‘) 

4-  C3k3sm3dcos0(l  4-  t'x2)4  -  C3k3sin0cos3<t>xA  (x2  _  2 z ) 2  } 
■  \-isin0(p4x8  4-  p3ix6  4-  P2X4  ~  Pi  tX2  4-  po)]- 
The  pressure  integral  along  the  SDP  becomes 

Psdp{R-0 )  =  V «)(pi  -  p2  -  p?J. 
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/  -= - r^dx  =  —erfc(bn)  ^  for  SR{6}  >  0. 
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Summing  Eqs.  (80)  and  (81)  for  the  k  —  0  terms  results  in 
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By  using  the  asymptotic  expansion  of  Hankel  function  given  in  Eq.  (66).  p0  is  then 


of  the  form 


=  £«  +  <) 


H()1\koRsin2<t>)( —  )erf  c(-ixn  \//knR)e  fc"Rv,‘ . 


For  integers  k  >  1.  substituting  Eqs.  (57).  (58)  and  (83)  into  Eq.  (80).  then  p]  can 
be  expressed  as: 
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To  obtain  an  exact  solution  for  an  observation  point  at  the  direction  normal  to 
the  plate,  one  may  have  to  start  with  the  integral  representation  of  the  pressure  field 


which  contains  the  Bessel  function.  Thus. 


Psdp{R,0)  = 
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where  u  =  koRsin0\cos<px V X2  -  2t  stnc>(l  4-  t x 2 )  1  - 

As  the  observation  angle  0  approaches  zero,  Ar(x,0)  and  /l,(x,<2>)  in  Eq.  (76) 
will  vanish,  and  Eq.  (89)  is  reduced  to 
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From  the  transformation  in  Eq.  (45).  it  can  be  seen  that  the  ten  simple  poles 
X^  are  now  reduced  to  five  poles  of  second  order  when  0  =  0°.  Thus,  the  partial 


fractions  for  the  integrand  becomes: 
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Substituting  Eq.  (91)  into  Eq.  (90)  and  changing  the  integration  variable  yields 
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The  remaining  two  integral  can  be  found  in  [16]  and  expressed  as 


--koRy 

/  - -dxj  =  - e-k"R Et{k0Rx2nh 

Jo  y-xl 


•OO  e~k„Ry 


Jo  = 

where  E,  is  the  exponential  integral  function  defined  as 
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Since  there  is  no  pole  located  between  the  oringinal  integration  path  and  the 
SDP  when  the  observation  angle  equals  zero,  the  total  pressure  can  be  written  as 


i  k  R  ^ 
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For  very  large  /cni?.  the  asymptotic  expansion  for  the  exponential  integral  function  is 


of  the  form 


Et{z )  ss  —  i  arg(-z)  i  <  3 7r / 2 . 


The  farfield  solution  for  the  radiated  pressure  at  the  direction  normal  to  the  plate 


becomes 
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where  a„  is  expressed  in  Eq.  (91). 


Notice  that  in  order  to  obtain  the  numerical  results  from  Eq.  (96).  it  is 
convenient  to  calculate  the  exponential  integral  function  by  using  the  Taylor  series 
for  small  arguments  and  the  asymptotic  series  tor  large  arguments.  The  Taylor 


series  is  given  as 


Ex(z)  =  -y  *  lo<j(-z)  - 


where  j  =  0.57721566-19  ••  -.  is  the  Euler  constant. 
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4.3  Numerical  Integration 

There  are  two  ways  to  perform  the  numerical  integration.  One  is  to  use  Eq. 
(71)  which  contains  the  double  integration  and  consequently  runs  considerable 
computer  CPU  time.  To  save  the  CPU  time.  Eq.  (71)  is  written  as 

p(R,<t>)  =  2^£2*L.(fcoW,-,/4, 

47rmn2 
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where  u  is  given  in  Eq.  (89).  The  Hankel  function  has  a  Taylor  series  expansion 


=  J0(z)+iY0(z), 


where 
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and  -v  is  the  Euler  constant,  and  the  asymptotic  expansion  for  large  arguments 
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It  should  be  noted  that  computing  the  Hankel  function  with  either  the  Taylor  or  the 
asymptotic  series  and  then  performing  the  single  integration  have  saved  the  CPU  run 
time  by  the  order  of  2  as  compared  with  those  for  the  double  integration. 


CHAPTER  5 


NUMERICAL  RESULTS  AND  DISCUSSIONS 


5.1  Introduction 

In  this  chapter,  numerical  results  are  presented  for  the  acoustic  radiated 
pressure  obtained  by  the  standard  steepest  descent  method,  the  modified  steepest 
descent  method  and  by  the  numerical  integration  for  both  line-  and  point-force 
excitations.  To  illustrate  the  influence  of  poles,  the  numerical  results  of  the 
asymptotic  series,  the  leading  term  approximation  or  the  numerical  integration  are 
compared  to  the  additional  contribution  due  to  the  poles  which  lie  between  the  SDP 
and  the  original  path  of  integration.  Comparison  is  also  made  for  plate  with  and 
without  structural  damping.  Some  additional  results,  such  as  behavior  of  integrand 
along  the  SDP.  surface  wave  pole  contribution,  grazing  and  nearfield  pressures, 
velocity  field,  and  acoustic  intensity,  are  also  presented. 

5.2  Pole  Contribution 

For  the  Timoshenko-Mindlin  plate  theory,  the  poles  can  be  found  by  solving  the 
following  equation  derived  from  the  denominator  of  the  integrand  in  Eq.  (15): 
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where  T2  =  K2  -  1  and  K  =  k.  kn. 

Extensive  discussions  on  the  poles  have  been  given  by  Stuart  [17J.  Crighton 
[18-19],  Pierucci  and  Graham  [20-21],  and  by  Strawderman.  Ko  and  Nuttall  [22] 
There  are  ten  simple  poles,  six  of  them  are  located  in  the  top  Riemann  sheet  and 
four  in  the  lower  Riemann  sheet.  The  typical  locations  of  poles  in  the  complex  I  - 
and  A'-plane  were  computed  numerically  and  are  shown  in  Fig,  5  The  locations  of 


poles  depend  on  the  frequency  only  and  Fig.  6  shows  the  frequency  loci  of  poles  in 
the  complex  A'-plane.  It  should  be  noted  that  if  the  plate  is  lossy,  i.e..  its  Young's 
modulus  is  complex-valued,  the  locations  of  poles  are  also  dependent  upon  the  loss 
factor  of  the  plate  material. 

There  are  only  three  poles  that  could  fall  between  the  SDP  and  the  original 
integration  path,  as  shown  in  Fig.  4.  Two  frequencies  are  of  importance,  which 
are  for  steel  plate  and  water  flc  =  1.23  and  flc0  =  6.79.  The  pole  03.  known  as 
the  leaky  wave  pole  due  to  the  flexural  waves,  is  included  when  fl  >  flc.  The  pole 
q4.  which  is  the  leaky  wave  pole  due  to  the  thickness  shear  mode,  is  included  when 
H  >  nco.  The  third  one  is  -ar>.  which  is  a  surface  wave  pole,  is  included  because  of 
the  addition  of  small  amount  of  damping. 

The  poles  affect  the  radiated  pressure  in  two  ways.  The  proximity  of  a  pole  to 
the  saddle  point  generates  a  peak  in  the  directivity  function  due  to  that  pole.  In  the 
SDP  solutions.  Eqs.  (41)  and  (68).  the  denominator  is  minimized  when  D  >  flc 
which  generates  a  peak  in  the  directivity  function  at  the  coincidence  angle  o  —  oc 
and  a  second  peak  at  o  =  oC0  when  the  frequency  U  >  nc0.  In  the  asymptotic 
series  solutions.  Eqs.  (59)  and  (86).  the  complementary  error  function  becomes  large 
when  the  argument  \n  is  of  small  value.  The  argument  \n  represents  the  distance 
of  the  pole  \n  from  the  origin  at  \  =  0.  which  represents  the  saddle  point.  Similar 
arguments  can  be  advanced  for  the  numerical  integration  scheme.  Eqs.  (50)  and 
(77-79),  where  the  integrand  becomes  large  as  the  variable  of  integration  \  comes 
close  to  one  of  the  poles  \n.  Thus,  all  three  solutions  generate  one  or  two  peaks 
in  the  directivity  function  when  one  or  two  poles  approach  the  saddle  point.  As 
mentioned  earlier,  the  pole  location  depends  on  the  frequency  only  and  the  saddle 
point  is  located  at  the  observation  angle  o.  Thus,  the  distance  between  the  pole  and 
the  saddle  point  depends  on  both  the  forcing  frequency  and  the  observation  angle 
To  examine  the  distance  between  the  poles  and  a.t  from  the  saddle  point. 


the  distance  |an  -  <p\  is  plotted  vs.  the  observer  angle  4>  for  O  =  1,5  and  10  in  Figs. 
7-9.  respectively. 

In  Fig.  7,  ft  =  1  is  below  the  first  coincidence  frequency  flc.  an's  distance 
is  large  for  all  observer  angle.  For  fi  =  5  >  f2c  in  Fig.  8,  Q3's  distance  becomes 
smallest  at  o  —  <t>c  =  39°.  while  a4  and  -Qo  s  distances  are  large  for  0°  <  d>  <  90°. 
At  the  frequency  H  =  10  >  fico.  Fig.  9.  aa's  distance  is  minimal  near  <pc  =  3-4° 
and  a4's  distance  is  smallest  near  qco  =  11°.  It  should  be  noted  that  the  higher  the 
frequency,  the  smaller  the  distances  of  03  and  a4  are  from  the  saddle  point.  This 
means  that  the  peaks  in  the  directivity  function  become  sharper  and  sharper  with 
increasing  frequency.  Since  the  poles  03  and  a4  have  limiting  values  when  Q  — *  00. 
the  coincidence  angles  have  high  frequency  limits  oc  —  30°  and  oco  — 1 •  17P  as 
f2  — »  oc  for  steel  plate  and  water. 

The  second  way  a  pole  influences  the  radiated  pressure  is  through  the  residues 
of  the  poles  located  between  the  SDP  and  the  original  integration  path.  These 
contributions  are  the  same  for  all  three  types  of  solutions.  As  can  be  seen  in  Eqs. 
(43)  and  (70).  an  exponential  decay  would  occur  for  a  complex-valued  a,,.  In  Eq. 
(43).  the  residue  decays  with  the  imaginary  part  of  cos(ari-(^).  It  can  also  be  shown 
that,  by  introducing  the  integral  representation  of  Hankel  function  expressed  in  Eq. 
(72).  the  residue  for  point  force  excitation  decays  with  /micos(an  -  <p)].  To  exhibit 
the  decay  rates  due  to  the  poles,  the  Imico.s(an  -  o) j  is  plotted  for  n  =  0,3  and 
•4  vs.  the  observer  angle  o  for  n  =  1.5  and  10  in  Figs.  10-12.  respectively.  If  the 
exponent  is  small,  then  the  residue  contributes  significantly  to  the  total  acoustic 
radiated  pressure.  For  H  <  6.79.  only  03  is  seen  to  have  a  small  exponent  near 
the  first  coincidence  angle  Or  while  the  exponents  for  the  other  two  poles,  a,,  and 
o4.  remain  relatively  large,  see  Figs  10-11  For  H  6.79.  a-  still  gives  a  small 
exponent  near  the  first  coincidence  angle,  but  now  a4  also  has  a  small  exponent 
near  the  second  coincidence  angle,  o  ,.  see  Fig  12.  Thus,  unless  the  multiplier  of 


Im\cos(an  -  <p)\,  koT.  is  very  large,  the  pole  may  contribute  significantly  to  the  total 
pressure,  especially  near  the  two  coincidence  angles  <pc  and  oco- 


5.3  Line  Force  Excitation 

5.3.1  Integral  Along  the  Steepest  Descent  Path 

The  integral  along  the  steepest  descent  path  is  given  in  Eq.  (49).  The 
integrand  has  an  exponential  term  e~kc,rx'  that  decays  faster  with  larger  value  of 
kor.  Thus,  the  integral  converges  faster  the  greater  the  observer  distance  from  the 
origin  and  the  integrand  resembles  a  Dirac-delta  function  at  the  origin  x  =  0  for 
increasing  kor.  This  indicates  that  the  major  contribution  to  the  integral  results 
from  the  region  near  the  origin  which  corresponds  to  the  saddle  point.  However, 
when  a  pole  is  located  near  the  saddle  point  or  on  the  SDP.  the  integrand  will  show 
another  peak.  Consequently,  the  presence  of  these  poles  requires  the  use  of  the 
modified  saddle  point  method  which  would  account  for  the  significant  contributions 
from  the  neighborhood  of  these  poles. 

Since  the  location  of  poles  in  relation  to  the  saddle  point  depends  on  the 
observation  angle  and  the  frequency,  one  expects  these  two  parameters  to  influence 
the  rate  of  decay  of  the  integrand  along  the  SDP.  For  frequencies  below  the 
coincidence,  e.g..  at  fi  =  1.0  and  d>  —  0°.  the  decay  rate  is  seen  to  depend  on 
the  value  kor  in  Fig.  13.  When  a  pole  comes  near  the  SDP.  such  as  the  case  when 
H  =  1  and  6  =  60°  in  Fig.  14.  one  must  account  for  the  contribution  of  the  pole 
to  the  integral.  At  the  grazing  angle  o  =  90°.  the  integrand  vanishes  at  the  saddle 
point.  Fig.  15.  However,  the  integral  itself  does  not  vanish  and  thus  requiring  the 
evaluation  by  an  asymptotic  series  instead  of  the  first  order  SDP  approximation.  The 
most  important  pole  is  the  leaky  wave  pole  which  comes  closest  to  the  saddle 
point  near  the  coincidence  angle  o-.  Thus,  the  integrand  becomes  more  like  a  Dirac- 
delta  function  at  the  coincidence  angle,  as  can  be  observed  at  o  =  o,-  =  31"  with 


n  =  10  in  Fig.  16.  When  the  thickness  shear  pole  q4  comes  close  to  the  saddle 
point,  a  similar  but  sharper  peak  occurs.  Fig.  17  for  o  =  oco  =  11°  and  Cl  =  10. 


5.3.2  Influence  of  Range  (kr>r)  on  the  Solution 

The  uniformly  asymptotic  series  solution.  Eq.  (59),  shows  the  dependence  of 
the  solution  in  terms  of  the  non-dimensional  parameter  k0r.  The  first  part  of  the 
solution  depends  upon  the  complementary  error  function  of  argument  -iXn\  k0r 
which  is  not  negligible  as  Xnv'^o'-  becomes  small.  However,  for  k^r  become  large, 
these  terms  generate  an  asymptotic  series  of  (kor)~(rn+1/2'> .  which  is  similar  to  the 
second  part  of  the  solution.  Thus,  when  kor  is  not  large,  the  dependence  on  k0r  is 
due  to  the  complementary  error  function  and  the  second  part's  1  j  y'kc,r.  However, 
for  kr,r  >  1  or  when  \n  is  not  small,  the  leading  term  of  both  parts  is  l/%//cor. 

5.3.3  Radiated  Acoustic  Pressure 

The  radiated  acoustic  pressure  was  computed  by  the  saddle  point  method 
(SPM).  Eq.  (42).  the  uniformly  asymptotic  series  solution  (UAS).  Eq.  (60).  and 
by  the  numerical  integration  (Nl).  Eq.  (47).  The  numerical  results  were  obtained 
for  frequencies  fl  =  1,5  and  10  and  k0r  ranges  from  10  to  10°.  In  this  section,  the 
numerical  results  and  discussions  will  be  focused  in  three  different  frequency  ranges. 

To  delete  the  range  dependence  of  the  radiated  acoustic  pressure,  the  directivity 
function  ( DF )  is  defined  as 

DF{dB)  =20 logip(r.o)  p( r.O0)'.  (104 

5. 3. 3.1  Frequency  Range  Below  First  Coincidence  H.-  =  1.23 

In  this  frequency  range,  the  radiated  acoustic  pressure  has  no  peak  in  the 
directivity  function,  since  the  poles  are  not  in  the  immediate  proximity  of  the  saddle 
point  Figs.  18a  and  19a  show  the  directivity  functions  for  ll  -  1  and  k..r  =  10  and 
100.  respectively.  The  solutions  obtained  from  the  three  methods  are  numerically 
very  close  over  most  angles,  and  the  difference  becomes  indistinguishable  for  knr  = 


1000  (not  shown).  The  contributions  of  the  poles  are  significant  only  when  k0r 
is  small  (nearfield).  Furthermore,  the  surface  wave  pole  — e*o  gives  non-negligible 
contributions  near  grazing  angles  o  =  90°.  Thus  the  surface  pressure  due  to 
the  solutions  without  the  residues.  Figs.  18b  and  19b.  the  SPM  solution  always 
vanishes,  while  the  UAS  and  Nl  solutions  predict  non-negligible  but  small-valued 
acoustic  pressure  at  o  =  90°.  As  k0r  increases,  the  surface  pressure  decreases 
even  for  the  UAS  and  Nl  solutions,  as  can  be  seen  from  Fig.  19b.  In  general,  the 
residue  contribution  to  the  total  pressure  seems  to  be  most  important  near  the 
surface  of  the  plate,  as  shown  in  Figs.  18c  and  19c  for  O  =  1  and  kor  —  10  and 
100.  respectively. 

5. 3. 3. 2  Frequency  Range  1.23  <  H  <  6.79 

For  frequencies  between  the  first  coincidence  flc  =  1.23  and  the  second 
coincidence  frequency  flco  =  6.79.  there  is  only  one  peak  occurring  at  <pc  in  the 
directivity  pattern.  This  peak  is  attributed  to  the  closeness  of  the  leaky  wave  pole 
c*3  to  the  saddle  point,  and  occurs  when  the  trace  acoustic  wavenumber  is  equal  to 
the  flexural  wavenumber.  Computations  were  performed  for  H  =  5  and  kor  =  10 
to  104  as  shown  in  Figs.  20-23.  For  kor  =  10.  100  and  H  =  5.  Fig.  20a  and  21a. 
the  approximate  SPM  solutions  predict  considerably  higher  acoustic  pressures  than 
the  more  accurate  UAS  and  Nl  solutions.  The  differences  in  the  directivity  functions 
between  the  predicted  SPM.  UAS  and  Nl  solutions  narrow  down  for  kor  =  103 
in  Fig.  22.  but  the  SPM  solution  is  markedly  different  from  the  other  two  near  the 
coincidence  angle.  The  differences  become  negligible  for  kor  =  104.  as  shown  in  Fig. 
23. 

The  role  played  by  the  residue  terms  (surrace  and  nearfield  effects)  are 
significant  when  one  compares  the  integral  along  the  SDP  in  Fig.  20a  with  the  total 
solution  in  Fig  20b  The  influence  of  the  poles  persists  even  for  k,,r  —  100.  shown 
in  Fig.  21a  and  21b.  and  produces  constructive  and  destructive  interference  near  the 


plate  surface.  The  contribution  of  poles  is  still  significant  over  the  angles  range  from 
40°  to  90°,  Fig.  22.  for  kor  =  103.  It  can  be  seen  that,  in  general,  the  higher  the 
frequency  and  the  lower  the  value  of  k0r.  the  more  influential  the  poles'  contributions 
become.  Figs.  20c  and  21c  show  the  poles’  contribution,  mainly  due  to  the  leaky 
wave  pole  0:3  for  0  =  5  and  k0r  =  10  and  100,  respectively. 

5. 3. 3. 3  Frequency  Range  Above  Second  Coincidence  f!f„  =  6.79 

For  frequencies  above  the  second  coincidence  frequency  Qco.  a  second  peak 
appears  in  the  directivity  function.  This  particular  peak  is  attributed  to  the  thickness 
shear  pole  0:4  being  dose  to  the  saddle  point,  i.e..  when  the  acoustic  wavenumber 
is  equal  to  the  structural  wavenumber  of  the  second  branch  of  the  plate  frequency 
spectrum. 


Computations  were  performed  for  fl  =  10  and  knr  =  10  to  10°.  For  low  ranges, 
knr  =  10  and  102  in  Figs.  24  and  25.  the  SPM  solution  generates  two  peaks  while 
the  UAS  and  Nl  solutions  have  none.  The  coincidence  peak  near  oc  =  34°  for  k0r  = 
103  in  Fig.  26  is  reasonably  modeled  by  the  SPM  solution,  but  the  second  peak  at 
0CO  =  11°  does  not  appear  yet.  Again,  the  pole  contributions  influence  the  directivity 
function  throughout  the  angular  range,  as  seen  from  the  constructive  maxima  and 
destructive  minima.  The  three  solutions  converge  at  kor  =  10s,  as  shown  in  Fig.  27. 
5. 3. 3. 4  Radiated  Pressure  at  <t>  —  0° 

In  the  previous  section,  the  directivity  function  is  computed  in  such  a  way  that 
the  pressure  at  any  aspect  angle  o  is  normalized  by  the  pressure  at  o  —  0°.  In  order 
that  one  can  obtain  the  absolute  value  of  the  radiated  pressure  due  to  the  line-force 
excitation,  the  radiated  pressure  at  o  -  0°.  normalized  to  .  is  plotted  in  Fig. 

28  for  fl  =  0.5.  1,  5  and  10.  Thus,  the  directivity  functions  in  Figs.  18-27  can  be 
corrected  by  the  addition  of  the  radiated  pressure  at  o  =  0"  shown  in  Fig.  28. 


The  radiated  acoustic  pressure  at  the  surface  of  the  plate  can  be  obtained  by 
setting  the  observer  angle  <p  equal  to  90°.  The  first  order  approximation  from  the 
standard  steepest  descent  method  vanishes  at  the  plate  surface,  as  can  be  seen 
in  Eq.  (41).  This  means  that  the  solution  depends  only  on  the  residues  of  the 
poles  located  between  the  SDP  and  the  original  integration  path.  Eq.  (43).  For  the 
asymptotic  series  solution,  the  leading  term  in  the  series,  m  =  0  in  Eq.  (59).  does 
not  vanish  in  general  and  contributes  a  term  of  the  order  l/y/ko?  to  the  radiated 
acoustic  pressure.  However,  the  leading  term  tends  to  be  zero  at  <j>  =  90°  as  k0r 
becomes  very  large.  This  means  that  the  remaining  terms  of  the  series,  m  >  1. 
contribute  a  term  of  the  order  ( kor)~ 3//2  for  kqr  >>  1.  Thus,  the  pole  contribution 
is  seen  to  be  very  important  to  the  solution  of  the  radiated  pressure  near  and  at 
grazing  angles. 

5. 3. 3. 5.1  Surface  Wave  Pole  and  Surface  Pressure 

For  an  observer  on  the  surface  of  the  plate,  i.e..  at  the  grazing  angle  <p  = 

90°.  the  surface  wave  pole  — a0  plays  an  important  role.  This  surface  wave  pole 
corresponds  to  -Kq  and  is  real  in  the  complex  wavenumber  ff-plane.  As  can  be 
seen  from  Fig.  6.  the  pole  —  Kq  is  always  greater  than  1  and  it  approaches  1  for 
H  —  oo.  This  means  that  -qq  =  rr / 2  -  «(fi).  where  e(Q)  >  0  is  a  real  positive 
number  that  depends  on  the  frequency  only.  The  pole  contribution,  expressed  in  the 
form  of  a  residue,  has  an  exponential  term  in  its  expression 

glfci»rC0J(  Q?1  —  <J>)  glknrcosi  z?  —  tp)coshl  ^  —  Knr3tn{  X 

_  zcosht  ^  —  k,,z si  nhi 

Thus,  when  the  observer  is  at  grazing,  i.e..  at  the  surface  of  the  plate  o  =  t  2 
or  ;  =  0.  the  exponential  term  in  Eq.  (105)  becomes  .  which  indicates 

an  unattenuated  wave  traveling  along  the  surface  of  the  plate.  Furthermore,  the 


—  -r>  1  nae  i 

(105) 


surface  wave  pole  -Qo  is  not  dependent  on  the  acoustic  medium  and  it  seems  to 
be  numerically  close  to  the  flexural  wave  propagating-mode  pole  of  a  plate  without 
acoustic  loading.  This  fact  can  be  gleaned  from  the  computation  of  the  surface  wave 
pole  for  both  air  and  water  loadings,  as  shown  in  Table  1.  It  can  be  seen  that  A'o(Q) 
or  e(f2)  does  not  differ  much  between  water  and  air  loadings,  despite  the  fact  that 
the  ratio  of  their  characteristic  impedances  is  as  high  as  3000. 

The  contribution  of  the  pole  -ao  in  Eq.  (105)  clearly  shows  a  surface  wave's 
behavior,  i.e.,  the  pressure  wave  propagates  unattenuatedly  along  the  plate  surface. 
z  =  0.  and  decays  exponentially  away  from  the  plate  surface  at  a  rate  of  k0zsinhe. 
Since  the  poles  03  and  a4  are  complex-valued  at  all  frequencies,  these  poles' 
contributions  to  the  surface  pressure  generate  exponentially  decaying  terms,  as  can 
be  seen  from  the  plots  of  /mlcosfa  -  <p)|  in  Fig.  10-12,  and  become  negligible  for 
large  k0r.  On  the  other  hand,  the  pole  — ao  generates  a  propagating,  unattenuated 
wave  traveling  along  the  plate  surface  at  all  frequencies  and  becomes  dominant  for 
k0r  large. 

To  investigate  the  role  played  by  the  poles'  contributions,  mostly  due  to  the 
surface  wave,  and  the  integral  contribution  to  the  total  pressure,  they  were  computed 
separately  and  plotted  for  H  =  0.5,  5,  and  10  in  Figs.  29 — 31 .  respectively.  In  these 
figures,  one  can  see  that  the  contribution  from  the  UAS  integral  in  Eq.  (59)  decays 
at  a  rate  (k0r)~a  on  the  surface,  where  kr>r  =  k0x.  The  decaying  rate,  which 
depends  on  the  frequency  and  k^r.  is  found  to  be  0.5  <  a  <  1.5.  and  the  upper 
limit  of  1.5  is  reached  for  k0r  >>  1  and  H  >>  1.  In  general,  since  the  integral  along 
the  SDP  and  the  contributions  of  poles  a3  ana  a4  decay  at  some  rate  with  k0r.  the 
total  field  would  also  decay  for  low  knr  and  then  become  constant  due  to  the  surface 
wave  at  some  distance  knr  away  from  the  line  force. 
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5. 3. 3. 5. 2  Nearfield  Pressure 

For  observer  located  away  from  the  surface  of  the  plate,  all  the  three  poles 
generate  exponentially  decaying  contributions,  as  can  be  seen  in  Eq.  (105).  The  rate 
of  exponential  decay  is  a  function  of  the  observation  angle  and  e.  and  the  latter  is 
listed  in  Table  1.  As  the  frequency  H  increases,  the  decay  rate  t  diminishes,  which 
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indicates  that  the  surface  wave  nearfield  extends  out  farther  and  farther.  However, 
since  ko  also  increase  with  fi.  the  decay  rate  kozsinhe  reaches  a  constant  limit  of 
0.05A:c2.  kc  being  the  acoustic  wavenumber  at  the  coincidence  frequency  nc. 
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As  mentioned  in  the  previous  section,  the  contribution  from  the  UAS  integral 
decays  at  a  rate  (A:or)-a.  ”a”  being  dependent  on  the  frequency  H  and  the  range 
kor.  thus  the  total  nearfield  pressure  would  be  dominated  by  the  UAS  integral.  For 
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low  frequency,  such  as  fl  =  0.5.  the  total  pressures  at  <p  =  87°.  88°  and  89°  are 
shown  in  Fig.  32.  It  can  be  seen  that  the  integral  (UAS)  becomes  more  dominant  at 


closer  ranges  k0r  with  decreasing  observer  angle  <p.  Similar  behaviors  were  exhibited 
at  higher  frequencies,  as  shown  in  Figs.  33  and  34  for  fi  =  5  and  10.  respectively. 


In  order  to  obtain  the  radiated  acoustic  pressure  near  the  surface  of  the  plate, 
the  radiated  pressure  at  a  distance  0.01A  away  from  the  plate  surface  is  plotted  vs. 
k0x  ranging  from  10  to  103  for  fi  =  1  in  Fig.  35.  Since  the  pressure  is  computed 
below  coincidence,  the  contributions  from  both  the  leaky  wave  pole  and  the  integral 
along  the  SDP  are  relatively  small.  Thus,  the  total  pressure  is  due  to  the  surface 
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wave.  As  the  observer  moves  away  from  the  surface  of  the  plate,  the  integral  along 
the  SDP  becomes  dominant  for  knx  increases  up  to  102.  as  shown  in  Fig.  36  for 
:  =  \.  For  higher  k0x.  the  integral  and  the  leaky  wave  decay  away  and  the  total 
pressure  could  be  attributed  to  the  surface  wave. 

For  frequencies  above  coincidence,  the  leaky  wave  pole  contributes  mostly  to 
the  radiated  pressure  for  knx  small,  while  the  integral  along  the  SDP  and  the  surface 
wave  are  insignificant,  as  shown  in  Fig.  37  for  H  =  5  and  c  =  0.01A.  However,  the 
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surface  wave  becomes  more  important  after  the  integral  and  the  leaky  wave  decay 
away  for  kox  >  103.  Similar  features  of  the  radiated  pressure  can  be  observed  in  Fig. 
38  for  ft  =  5  and  z  =  A. 

5. 3. 3. 6  Influence  of  Structural  Damping 

The  Young's  modulus  encountered  in  the  computation  so  far  is  considered 
real-valued,  i.e..  the  plate  material  is  purely  elastic.  In  this  section,  the  influence  of 
the  added  structural  damping  on  the  locations  of  poles  and  the  radiated  acoustic 
pressure  will  be  discussed.  The  absorption  of  the  elastic  energy  can  be  taken  into 
account  by  using  the  complex  Young  s  modulus  which  is  related  to  the  real  one  by 
the  following  equation 

E  =  E(  1  —  iT7 ) ,  (106) 

where  r/  is  the  plate  material  loss  factor  and  0.0  <  rj  <  1.0  for  nearly  all  materials  of 
interest.  For  small  amount  of  damping,  the  sound  speed  is  corrected  by  multiplying 
the  sound  speed  by  (1  -  trj)l/ 2  as  (1  -  ^ir?).  and  hence  the  coincidence  frequency  in 
Eq.  (17)  becomes  wc(l  4-  ^ir 7). 

Taking  into  account  the  dissipation  of  elastic  energy  by  introducing  a  complex¬ 
valued  Young's  modulus  usually  causes  all  the  poles  an  to  rotate  clockwise  in 
the  complex  wavenumber  plane.  It  should  be  noted  that  the  leaky  wave  poles  are 
defined  as  those  located  on  the  bottom  Riemann  sheet  of  the  complex  K- plane  and 
must  appear  near  the  saddle  point  at  some  frequencies.  For  given  plate  properties 
and  dimensions,  the  leaky  wave  pole's  distance  from  the  real  axis  is  a  function  of 
frequency  and  is  dependent  on  the  plate's  absorption  of  the  vibrational  energy,  i.e.. 
on  the  magnitude  of  the  loss  factor. 

Below  coincidence,  since  the  energy  radiated  into  the  surrounding  medium  is  due 
only  to  the  line  force,  the  radiated  acoustic  pressure  shown  in  Fig  39  for  k,,r  =  10 
and  H  =  0.5.  displays  no  significant  dependence  on  the  structural  damping  in  the 


whole  fieid.  In  the  high-frequency  range  above  coincidence,  the  radiation  pattern  is 
significantly  affected  by  structural  damping  because,  in  addition  to  the  external  force, 
the  whole  plate  also  radiates  energy.  As  can  be  seen  in  Fig.  40a  for  fi  =  5.0  and 
kor  =  102.  the  nearfield  pressure  attributed  to  the  leaky  wave  pole  is  remarkedly 
reduced,  while  the  influence  on  the  integral  along  the  SDP  near  the  coincidence  angle 
is  relatively  small,  as  shown  in  Fig.  40b.  The  influence  of  structural  damping  on 
residue  contributions  and  the  integral  along  the  SDP  can  be  explained  from  Eqs. 

(42)  and  (59).  For  small  kor.  the  u>-function  in  Eq.  (59)  is  evaluated  by  the  power 
series  expansion,  Eq.  (61).  and  has  the  approximate  numerical  value  of  1.  Thus,  the 
contribution  due  to  the  integral  is  mostly  proportional  to  1  /\n.  Consequently,  the 
influence  of  structural  damping  is  not  significant  because  of  the  slight  change  in  \n- 
On  the  other  hand,  slight  change  of  the  pole's  distance  from  the  saddle  point  would 
result  in  a  faster-decaying  wave  since  the  residue.  Eq.  (43).  contains  an  exponential 
term  e_fcnrx" . 

For  farfield  pressure.  k0r  =  10°  at  fl  =  5.  the  contribution  from  the  leaky 
wave  pole  becomes  insignificant  and  the  only  influence  of  damping  on  the  integral 
along  SDP  appears  near  the  coincidence  angle,  as  shown  in  Fig.  41.  The  radiated 
pressure  is  significantly  reduced  due  to  the  fact  that  certain  amounts  of  the  plate's 
vibrational  energy  is  dissipated  through  the  damping  mechanism.  Similar  phenomina 
can  be  observed  at  higher  frequencies  in  Figs.  42  and  43  for  fl  =  10  and  k0r  —  102 
and  10°.  respectively. 

At  grazing  angle,  the  surface  wave  pole.  Kn.  is  no  longer  located  on  the  real 
axis  and  shifted  away  from  the  real  axis  in  the  /\ -wavenumber  plane  because  of  the 
added  damping.  Thus,  the  surface  wave  generates  3  decaying  term,  see  Figs.  44  and 
45.  Figs.  46-47  show  the  influences  of  different  amount  of  damping  on  the  grazing 
field  for  H  =  0.5  and  5.0.  respectively. 


5.4  Point  Force  Excitation 


5.4.1  Integral  Along  the  Steeoest  Descent  Path 


For  point  force  excitation,  the  pressure  integral  along  the  SDP  in  given  in  Eq. 
(100).  The  integrand  also  contains  an  exponential-decaying  term  e~k"Rx~ .  Thus, 
as  k0R  increases,  the  integral  converges  faster  and  resembles  a  Dirac-delta  function 
at  the  origin  x  =  0  which  corresponds  to  the  saddle  point.  Since  the  numerator 
of  the  integrand  consists  of  both  even  and  odd  functions  in  x.  and  only  the  even 
part  will  contribute  to  the  pressure  integral,  the  integrand  is  evaluated  as  one  half 
of  the  sum  of  the  integrand  with  both  x  and  — x.  As  can  be  seen  in  Figs.  48-52.  the 
same  features  as  for  the  line  force  excitation  are  observed  when  the  leaky  wave  pole 
appears  near  the  saddle  point.  The  notches  that  appear  in  these  figures  are  due  to 
the  Hankel  function  s  local  minima. 

5.4.2  Radiated  Acoustic  Pressure 

The  expression  for  radiated  acoustic  pressure  for  point  force  excitation.  Eq. 

(86).  is  similar  to  Eq.  (59).  which  is  for  line  force  excitation.  Thus,  most  of 
discussions  which  are  same  for  line  force  excitation  will  not  be  given  in  this  section, 
but  can  be  referred  in  detail  in  chapter  5.3.  Since  the  results  have  shown  that  the 
directivity  functions  are  exactly  same  as  those  for  line  force  excitation,  thus  only  the 
pressure  at  <p  =  0°  is  needed  and  shown  in  Fig.  53.  The  directivity  functions  shown 
in  Figs.  18-27  can  then  be  corrected  by  Fig.  53  to  obtain  the  radiated  acoustic 
pressure.  It  is  observed  that  the  decay  rate  for  the  radiated  pressure  at  o  =  0°  is 
(kn  R)~2. 

5.4.3  Grazing  and  Nearfield  Radiated  Pressure 

At  the  surface  of  the  plate,  o  =  90°  and  k^R  =  k,,r.  the  residue  of  the  surface 
wave  pole  -an  =  ~  2  -  t'e(fl)  contains  a  Hankel  function  as  follows: 


//0*  ,knRsin(  —  «o)i  —  Hlnl){knRcosht), 


(106) 


/..VjV1 


•  WV 

•vv\ 

“WV 

iasii 

vw 

•  ■ 

*  /  v  Vv 
»  V  V  V 

»  rv  v 
t  j*  , 


.  -v«v 

.vv\ 

',vvv: 


m 


V  A 


-Vv' 

' A 
A/  v  ■ 


/  V  V 


& 


*  ■ » ^  * 1 

■  vV'.' 

■  V^v-V 

y.v.vv 

\  f««.  A 


-vv.y, 
■>  v 

M' 


% 

•V%V 
•*>>> 


tt  vr  Tx  V  V".  . 


a 


where  £  >  0  is  dependent  on  the  frequency  only.  The  poles'contributions  and  the 
integral  contribution  to  the  total  pressure  at  the  grazing  angle  are  shown  in  Figs 
54-56  for  H  =  0.5.  5  and  10.  respectively.  As  can  be  seen  from  these  figures,  the 
contribution  from  the  integral  along  the  SDP  decays  at  a  rate  lk,iR)  1  and  "a"  is 
found  to  be  1.0  <  a  <  2.0,  and  the  upper  limit  is  reached  for  k,,R  1  and 

H  >>  1.  It  should  be  noted  that  the  radiated  pressure  due  to  the  surface  wave  is 
no  longer  constant  near  the  surface  of  the  plate  as  was  seen  in  Figs  29-31  for  line 
force  excitation,  but  decays  at  a  rate  (k0R)~l/2.  Figs.  57-62  show  the  nearfield 
radiated  pressures  for  fl  =  1.  5.  10  and  z  —  0.01A  and  A,  respectively.  Similar 
discussions,  except  for  the  decay  rates  of  UAS  integral  and  the  surface  wave,  have 
been  given  in  chapter  5.3  for  line  force  excitation. 


CHAPTER  6 


CONCLUDING  REMARKS 


6.1  Summary  and  Conclusion 

A  new  solution  to  the  acoustic  radiation  from  an  infinite  fluid-loaded,  elastic 
plate  excited  by  a  line  or  a  point  force  has  been  obtained,  for  the  first  time,  in  the 
form  of  a  uniformly  asymptotic  series.  It  was  shown  that  the  uniformly  asymptotic 
series  solutions,  derived  from  the  modified  steepest  descent  method,  indicate  clearly 
that  the  first  order  approximation  obtained  previously  by  other  researchers  using  the 
standard  saddle  point  method  overpredicts  the  acoustic  pressure  in  the  nearfield. 
Furthermore,  it  was  shown  that  the  error  in  predicting  the  acoustic  pressure  using 
the  saddle  point  method  is  especially  pronounced  at  the  coincidence  and  grazing 
angles  when  compared  to  that  due  to  the  uniformly  asymptotic  series. 

It  was  also  shown  that,  rather  than  performing  the  original  integral  with  branch 
cuts,  the  transformed  integration  along  the  steepest  descent  path  (SDP)  is  more 
efficient  numerically.  The  excellent  agreement  between  the  two  solutions,  the 
uniformly  asymptotic  series  solution  and  the  numerical  integration  performed  along 
the  SDP.  validates  the  use  of  the  modified  saddle  point  method  for  predicting  the 
acoustic  pressure  in  the  nearfield  and  near  the  grazing  angle.  However,  the  widely 
used  standard  saddle  point  method  agrees  with  the  other  two  solutions  only  in  the 
distant  farfield  (e.g..  kor  >  104  for  fi  =  5.0). 

The  radiated  nearfield  acoustic  pressure  was  accurately  derived  in  this  study  and 
was  shown  to  be  dominated  by  the  UAS  (integral  along  the  SDP)  only,  which  decays 
at  the  rate  ( kor)~a .  where  "a"  ranges  from  0.5  to  1.5  for  line  force  excitation, 
depending  on  the  range  kor  and  the  frequency  H.  The  surface  pressure  was  shown 


to  be  dominated  by  the  surface  wave  poie.  which  represents  a  free  flexural  wave 
propagating  along  the  surface  of  the  plate  without  fluid  loading. 

The  influence  of  the  structural  damping  on  the  acoustic  radiation  has  also  been 
investigated.  In  general,  below  coincidence,  where  the  radiated  pressure  is  primarily 
due  to  the  nearfield  of  the  external  force,  the  added  damping  has  negligible  effect. 

For  frequencies  above  coincidence,  the  radiated  pressure  is  significantly  affected  by 
the  structural  damping  since  the  plate  as  well  as  the  region  near  the  external  force 
radiate  energy.  In  the  nearfield  and  above  the  first  coincidence  frequency,  where  the 
leaky  wave  dominates,  the  inclusion  of  structural  damping  shifts  the  leaky  wave  pole 
farther  away  from  the  saddle  point  and  hence  generates  more  exponential  decay. 
However,  the  integral  along  the  SDP.  which  depends  linearly  on  the  pole  location, 
is  slightly  reduced  only  near  the  coincidence  angle.  In  the  farfield.  where  the  solution 
due  to  the  leaky  wave  pole  decays  rapidly  and  the  integral  along  the  SDP  becomes 
dominant,  the  effect  of  damping  is  significant  only  near  the  coincidence  angle.  The 
structural  damping  also  shifts  the  real  surface  wave  pole  away  from  the  real-axis 
in  the  complex  wavenumber  plane  and  hence  generates  a  decaying  surface  acoustic 
pressure  which  depends  on  the  amount  of  damping  added. 

In  conclusion,  the  new  solutions  derived  in  this  thesis  in  the  form  of  a  uniformly 
asymptotic  series  and  a  fast-convergent  numerical  integral  were  shown  to  predict 
the  radiated  acoustic  pressure  accurately  in  the  nearfield  and  at  the  grazing  angle  in 
contrast  to  solutions  obtained  previously  by  first  order  approximations.  Since  the 
difference  in  the  predictions  persists  for  non-dimensional  distance  fcnr  over  10'5  to 
10T  the  uniformly  asymptotic  series  solution  cotained  in  this  thesis  should  be  used 
m  many  instances  of  practical  interest. 

6  2  Future  Work 

The  analysis  in  this  thesis  has  been  concentrated  on  the  acoustic  radiation  from 
infinite  homogeneous,  isotropic  elastic  plates  excitea  by  a  line  and  a  point  forces. 


However,  with  the  application  of  the  Green's  function  theory,  the  analysis  can  be 
easily  extended  to  handle  the  same  problem  but  with  arbitrary  disturbance.  It  can 
also  be  applied  to  the  nearfield  acoustic  radiation  from  a  rnoment-excited  plate. 
Further  work  of  interest  remaining  to  be  done  are  the  analyses  for  finite  plates  with 
boundary  effects,  as  well  as  non-isotropic  or  inhomogeneous  plates. 


Line  Force  F{t)6(x) 


Vacuum 


Figure 


or 

Point  Force  F(f )  ~~ 


Geometry  of  a  line-  or  point-excited  fluid-loaded  elastic 
plate  and  the  coordinate  system. 
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-a,,;  true  surface  wave  poles. 

— a i.4 :  complex  surface  wave  poles, 

=ar3:  leaky  wave  poles(  also  ran  rj  Q  >  fl,„  ). 


igure  3  Typical  locations  of  poles  in  the  complex  a-plane. 


Bottom  Riemann  Sheet 


Top  Riemann  Sheet 


Figure  5  Typical  locations  of  poles  in  the  complex  Y-  and  A'-plane 
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Figure  7  Distance  between  the  saddle  point  and  the  poles  in  the 
complex  a-plane.  H  =  1.0. 
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Figure  10  Absolute  values  of  Im  cos{a  -  o)\  H  =  1.0. 
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Figure  11  Absolute  values  of  Im  cos(a  -  o)'.  fl  =  5.0. 
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X.  distance  from  saddle  point 

16  Magnitude  of  integrand  along  SDP-  line  excited  plat 
at  0  10.0  and  4>  3-I1’. 
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Figure  20  Continued 
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Magnitude  of  integrand  along  SDP-  point-excited  plate 
at  H  =  100  and  4>  =  31°. 
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Radiated  pressure  from  a  point-excited  steel  plate 
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APPENDIX 


THE  TAYLOR'S  SERIES  EXPANSION 


In  Eq.  (77).  the  second  integral,  but  not  the  exponential  term,  can  be  expanded 
into  Taylor  s  series  in  x  and  y  about  (x,y)  =  (0,0).  Since  the  integrands  with  odd 
power  term  in  x  will  not  contribute  to  the  radiated  acoustic  pressure,  only  those  with 


even  power  in  x  are  considered.  Thus. 


y(x,y)\/l  +  T 
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where  y(x,y)  is  given  in  Eq.  (74)  and  fki  is  the  Taylor's  series  coefficient.  For 
reference,  the  coefficients  for  an  expansion  through  x4y4  are  given  as  follows: 
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Similarly,  the  second  integrals  in  Eqs.  (78-79)  are  expanded  3S 


i  -  1X1 
_ 2_ 

y(x>y) 


—  V 2e  ' 4  V]  9kix2ky‘l- 

k =0  1=0 


A  -  3) 


*  ’•  /'  «  V«  ’  7  '  •  V  ’  V  *  ■  V  ' 


v'VVV 


*  j 

*.v 

>KV 

\  1 


v- 

>.  v 
>^.v> 

.AV.V.V 

vx-M 


v‘.Nr.N'^ 


•  .N  . 


*  V  *. 

.*  V  V  V 


V  V  V 

•  '  »  Vv  »  v 
V  \s'  v*^ 

-V-\  .\V 

■  V  V  *JI  fyl 


#  v-  v  v 


A 


.1  i.t 


64  s:n°<£ 
i(l  +  2cos24>) 


5  +  10co52(6 


t(27  4-  78cos2©) 


256 sin' 0 
7  4-  28 cos26 
64  sin5<p 

—  t(47sin4<6  4-  330 sin2  0cos2  0  4-  315cos4<j>) 

256sin~<p 

—  (357sm4<z>  4-  3l50stn2<j>cos2<£  4-  3465cos4<i>) 

2048 sing0 

OO  OO 

-r^^EtaV. 

»<*■»>  fes 


(yi  -  5) 


-icosd) 

10  — 

2sin26  ’ 
— 3cOS0 

11  = 

8 sm40  ’ 
t locos® 

12  = 

64sine0 

—bcosd 

20  = 

8sm*o  ' 

(A -6) 


32sznv  o 

co.so(195  —  120cos2o) 


obsinBo 


